Draft version August 13, 2012 

Preprint typeset using I^-T^X style emulatcapj v. 5/2/11 



(N 

o 

< 
O 



X'- 

a. 



(N 
> 
O 

ov 

06 
o 

(N 



13 



THE EXTREMELY HIGH PEAK ENERGY OF GRB 110721A IN THE CONTEXT OF A DISSIPATIVE 
PHOTOSPHERE SYNCHROTRON EMISSION MODEL 

Peter Veres 1,2,3 '*, Bin-Bin Zhang 1 , Peter Meszaros 1,2,3 

(Dated: August 13, 2012) 

Draft version August 13, 2012 

ABSTRACT 

The Fermi observations of GRB 110721A (jAxelsson et al.l [2012T l have revealed an unusually high 
peak energy ~ 15 MeV in the first time bin of the prompt emission. We find that an interpretation is 
unlikely in terms of internal shock models, and confirm that a standard black-body photospheric model 
also falls short. We show that dissipative photospheric synchrotron models ranging from extreme 
magnetically dominated to baryon dominated dynamics, on the other hand, are able to accommodate 
such high peak values. 



1. INTRODUCTION 

High energy (> 100 MeV) observations of gamma-ray 
bursts (GRBs) have gained renewed inte rest since the 
launch of Fermi LAT ([Atwood etaD 120091 ) . This instru- 
ment has revealed a curious diversity in the high-energy 
behavior of GRBs and c onsequently repr esents a new 
challenge for models (see [Meszaros 201 2], for a r ecent 
review). For GRB 110726A (jAxelsson etal.1 [20121 ) . the 
extremely high e pca k ~ 15 MeV peak energy (peak of the 
eF e spectrum) , which was observed right after the onset 
of the burst, makes its interpretation in the framework of 
a simple internal shock synchrotron model challenging. 
Interpreting the peak as a modified black-body from a 
simple photosphere is also difficult, as it can only account 
for energies up to few MeVs, and even the subsequent 
time-bins show peak fluxes which, even though lower, 
are still unusually high. Such black-body photo sphere 
models have been considered bv lFan et al.l ( 20121 ) to ex- 
plain a L oc £p C 4 ak correlation found by iGhirlanda et al.l 
(|2012t ). but as shown by I Zhang et~aT1 (|2012D . they must 
lie below a line in the £ pca k — L plane which excludes the 
observed values for GRB 110721 A. 

Here we show that the high peak energy and lumi- 
nosity of GRB110721A can be interpreted in the frame- 
work of a different class of photospheres, where dissipa- 
tion occurs near the photospheric radius, and the spec- 
trum is characterized by a synchrotron peak. Dissipa- 
tive synchrotron photospheres have b een considered for 
baryo ni c dominated regimes by, e.g. iRees k, Meszaros! 
((20051) : IPe'er et al.l (1200(1 : IBeloborodovl (120101) and 1m 
magnetically domin a ted r e gimes by IGianniosI ([2001 : 
iTchekhovskov et all ([20121 ); IGianniosI (|2012}T Here we 
consider such photospheres for an arbitrary acceleration 
law, characteristic of either a baryonic or a magnetic 
dominated regime. Also the observed thermal bump 
close to 100 kcV is naturally incorporated in this model. 



1 Department of Astronomy and Astrophysics, Pennsylvania 
State University, 525 Davey Lab, University Park, PA 16802, 
USA 

2 Department of Physics, Pennsylvania State University, Uni- 
versity Park, PA 16802, USA 

3 Center for Particle Astrophysics. Pennsylvania State Univer- 
sity, University Park, PA 16802, USA 

Email: vcrcsp@psu.edu 



2. DYNAMICS, NON-DISSIPATIVE PHOTOSPHERES AND 
INTERNAL SHOCKS 

The initial acceleration behavior of the jet material is 
taken to be given by 



r(r) oc 



const. 



if r < r sat 
if r sat < r, 



(1) 



where r sat is the saturation radius beyond which the flow 
reaches its asymptotic coasting value rj = (L) / (Mc 2 ) = 
Pfinai- Based on the a nalysis of (e.g. IDrenkhahnl 120021 : 
lMeszaros~fc Recs 2011), for simple conical magnetically 
dominated models \x = 1/3, while in the simple baryoni- 
cally dominated case /i = 1 . More general cases involving 
different magnetic geometries or flows where the opening 
angle varies with radius will generally lie between 1/3 
an d 1. The case of \i = 1/3 has been further developed 
by iVeres fc Meszarosl (|2012t ) , including the effects of the 
shocked electrons at the external deceleration radius; the 
generalized magnetic model including an arbitrary accel- 
eration law 1/3 < // < 1 is considered in more detail in 
Veres et al. 2012 (in preparation). Here we calculate 
the relation between the luminosity and the peak energy 
based on the dissipative synchrotron photospheric model 
for values of between the two extremes, in the context 
of GRB 110721A. We denote physical quantities through 

the usual Q x = Q /10 x notation . 

As discussed by I Zhang et al.l (|2012f ) in the usual \i = 1 
baryonic dynamics the photospheres usually arise in the 
r > ^sat region, and in the absence of dissipation the 
black-body peak energy and flux representing a putative 
Band spectrum peak would fall below a line which ex- 
cludes the observations of GRB 110721A. 

Internal shocks may also be considered as a possible ex- 
planation, since this model is usually used for interpret- 
ing t he prompt gamma-ray emission ([Rees fc Meszarosl 
1!)!) 1 ;. For variability timescalcs t v seconds and an 
average Lorentz factor rj the internal shock occurs at 



a radius ns 



rfct v 



5 x 10 £^-1772 5 cm 



and for the usual magnetic field prescription one ex- 
pects an observer frame peak energy of ep^ak ~ 

1 ^52 2 ^5^-l e B,U.3 e e,0.5T|i McV ' &U U PP Cr limit 

based on a variability t v = 10" 1 s (the light 
curve of GRB 110721A is rather smooth, FRED- 
likc, suggesting a longer variability time). A some- 
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what different cal culation including pair formation 
(|Guetta et al.l 120011 ) indicates an upper limit e pca k < 

1 L^^^B^^vt^^v-aTTz MeV which is more strin- 
gent. Thus, standard internal shocks cannot explain this 
burst. 

3. A DISSIPATIVE PHOTOSPHERE MODEL 

Photospheric radius We consider a dissipative photo- 
sphere, where a substantial fraction ~ 0.5 of the bulk 
kinetic energy is converted into random particle energy. 
An extra spectral component coming from the interac- 
tion of the photospheric photons with shock-accelerated 
electrons at the external deceleration radius is neglected 
here, as this produces a GeV range contribution, whereas 
110721A has only a tentative GeV component. The pho- 
tosphere arises at optical depth r = 1, in a jet which 
initially accelerates as: T = (r/ro)**. Here ro is launch- 
ing radius at the base of the jet. The value of the pho- 
tospheric radiu s can be expressed in te rms of a critical 
Lorentz factor (|Meszaros fc Rees 2000) which, general- 
ized to an arbitrary is given by 



(La T /8irm p c 3 r Q )»/ 1+:if *. 



(2) 



For rj > r)T, the photosphere occurs in the accelerating 
phase, r p h = r^^{rjT /ri) x ^ 1+2ii \ while for r\ < rjr, it 
is r ph — Vt lfJ (jlr I v) 3 ■ The saturation radius is r sat = 

The photosphere occurs in the acceleration phase if 
T] > r/T, which is typical for a magnetically dominated 



(/i = 1/3) case, where t\t — 120 L 



1/6-1/6 

53 '0,7 



On the 



other hand, the photosphere is in the coasting phase for 
rj < tjt, which is typical for baryonic cases (fi = 1), where 

r/T — 1300 L^r Q The photospheric radius can be 
increased by a factor of ~ few by the presence o f pairs 
(|Bosniak fc Kumar! [20lll Veres fc Meszarosll2012l ). 

Synchrotron peak We expect that about half of the 
total luminosity wil l be converted to radiation close 
to the photosphere (|Gianniosll2012| ). The peak energy 
will form at moderate optical depths r > 1. The 
dissipation is expected to occur in mildly relativistic 
shocks with relative Lorentz factors of r r > 1, the ran- 
dom Lorentz factor of the shock-accelerated electrons is 
7o, P h ~ e e (m p /m e )T r « 600 e e , .5r r ,o, where e e is the 
fraction of the energy in electrons. The peak synchrotron 

energy can be expressed as e poak = i|^7o,phTf¥' 
where -B' p h = (327resWpC 2 n'b) 1 / 2 r r . The physical con- 
stants have the usual meaning, < 1 is the fraction of 
the energy in magnetic form, n' b = L/47rrp h m J ,c 3 F p h?7 is 
the comoving baryon density. This latter quantity scales 



as r 



up to the saturation radius, and as r in the 
coasting phase. The Lorentz factor of the photosphere is 
r p h = (fphAoy if the photosphere occurs in the acceler- 
ation phase, and it is n otherwise. All this put together 
results in a peak energy dependence 

£peak oc | L^rj~^r^e 2 e Tf/(l + z) if V > rjT ( 3 ) 

pcak \i-vv^r»/(i + z) if v < VT . 

We provide exact values for the representative cases dis- 
cussed here in Section [4] 



Thermal peak Besides the synchrotron component we 
expect also a (modified) thermal component from the 
photosphere. This component is observed as a ther- 
mal peak and it is advected from the launching radius, 
and cooled according to the generalized dynamics. We 
present concrete values for the temperature in Section 21 
here we only show the general dependence on the param- 
eters: 



T obs (r ph )cx { n >o 



£-5/12^8/3^/6/^ 



/{I + z)if T) > 77^4) 
if 7? < r/T- 

There were hints of a separate thermal component below 
the Band p ea k dPage et al.l I20TH: iGuiriec et al.l 120111: 
Zhang et afl |201lt IMcGlvnn fc Fermi Collaborationl 
2012D in previous bursts, and there are reportedly 
significant thermal c omponents also in GRB 110721A 
(|Axelsson et al.ll2012D . 

4. EXTREME MODELS 

Here we consider the £ pea k — L pairs for different \x 
values to show that in the framework of a dissipative 
photosphere, the high peak energy of GRB110721A can 
be obtained in a straightforward manner. 

From Equation [3J r\ > tjt case (photosphere in the ac- 
celeration phase), one sees that higher luminosities, and 
lower Lorentz factors and launching radii, will result in 
larger peak energies. Increasing L and decreasing ro will 
increase rjr eventually leading to the 77 < rjr case (pho- 
tosphere in the coasting phase). This transition is shown 
by the break in the evolution of e pca k with luminosity (see 
Figures [T] and [2] and note that in the extreme magnetic 
case, the peak energy does not depend on the luminos- 

ity). 

Figures [T] and [2] show that the observed highest peak 
energy- luminosity pair for GRB 110721A is in the ad- 
missible region of the diagrams, and for other reasonable 
parameters the peak energy is comfortably within the 
operating range of these models (admissible regions are 
under the modeled broken power laws, represented by 
dotted, dashed and dash-dotted lines). 

Here we present cases to illustrate the peak energy de- 
pendence of the luminosity, and show that within dissipa- 
tive synchrotron photospheric models we can reproduce 
the high peak energy. We have taken a redshift z = 0.382 
(|Bergerlf2011h for the calculations presented here. The 
luminosity values in the time bins are scattered around 
L = 10 52 erg/s (see e.g. Figure [T]) thus we take this as a 
reference value. 

4.1. Baryonic photosphere model (fi = 1) 

In this model, the magnetic field is subdominant, the 
dynamics are governed by the baryons in the outflow. 

The critical Lorentz factor is ryr ~ 740 L\j^r \^ , which 
puts the photosphere in the coasting phase for moder- 
ately high 77 < 600 values. The peak energy will become: 



e poak w 16 W^Ll^o.!^-! MeV. This is the right 
order of magnitude for not too extreme parameter val- 
ues. Lending further support for this model (the bary- 
onic variant) is the temperature of the thermal compo- 
nent, which also turns out the right order of magnitude, 

r« 70 L-^Mr keV. 



,1/2 
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4.2. Extreme magnetic model (/j, = 1/3) 

In this extreme case the break energy does not de- 
pend on the luminosity, nor the coasting Lorentz fac- 
tor. Still, for reasonable values, we can get (see equa- 
tion [3]) a peak value £ pca k ~ 13 r^J ^r.o.5 e B o 

McV. 

The temperature of the blackbody radiation is T ps 
3 -£'52 1 ^ 60r ?2. / 5 L5r o~7^ 30 k°V. The dependence on the pa- 
rameters is weak and this model cannot account for the 
observed temperature of the black-body component. 

4.3. Moderate magnetic model (fx = 0.5) 

Here magnetic fields are dominant, but the role of 
baryons is more marked than in the extreme magnetic 
case. The peak energy in this intermediate variant of 

the model is e pC ak ~ 14 £52%2^ /8r ^7 /8r ?,u.i4(o MeV. 
Again, this model can incorporate the rather large peak 
energy for not too exceptional parameters. The tempera- 
ture of the black-body is lower than in the baryonic pho- 
tosphere case, T sa 24 lI^Vx^o \^ kcV, but might 
still be consistent with the observations. 




1Q 48 1Q 50 1Q 52 1Q 54 

Luminosity [erg/s] 



Fig. 1. — GRB 110721A values of spectral peak energy versus 
peak luminosity for di fferent time bins ad vancing from top to bot- 
tom (diamonds, from [Zhang et al. (2015)). compared to expected 
limiting values from dissipative synchrotron photosphere models 
(dashed: baryonic fj, = 1, dotted: extreme magnetic = 1/3, dot- 
dashed: moderate magnetic fi = 0.5 (this work)). These lines are 
an upper limit: values of £ pea k an d L below them are allowed. 
Thus, moderate magnetic and baryonic dissipative photospheres 
are compatible with the highest peak energy values observed, and 
even the extreme magnetic case is within a standard deviation. 
The parameters used are: for fj, = 1/3: ro = 10 7 cm, eg = 1, 
T r = 2, for fi = 0.5: r = 10 7 cm, e B = 1, T r = 1.2, for fi = 1: 
r = 10 8 cm, e B = 0.1, T r = 1.2, and r/ = 300 throughout. The 
solid g ray line is a stand ard (non-dissipativc) photosph ere curve fol- 
lowing H*'an et af] 120121) . confirming the conclusion of lZhang et al.l 
120121 ) that non-dissipative photospheres cannot explain the high- 
est peak energies. 
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Fig. 2. — Same notation and parameters as Figure[T] except for 
the coasting Lorentz factor, which is r\ = 600 here. 

of iZhang et al.l <|2012f ) that a (non-dissipative) standard 
black-body baryonic photosphere model also cannot ex- 
plain such high peak values and fluxes. However, we show 
that dissipative photosphcric models, with a typical peak 
energy due to synchrotron radiation, arc able to accom- 
modate such high peak energies and flux values, with 
reasonable parameters, for cases where the dynamics is 
either baryonically or magnetically dominated. If the 
temperature of the putatively observed black-body com- 
ponent can be used as a discriminant, this would seem to 
favor a more baryon dominated dissipative photosphere 
model, fi closer to 1, although a moderate magnetically 
dominated photosphere may also be possible. 
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5. DISCUSSION 

We have considered the high peak energy values of 
the prompt spectrum of GRB 110721A, which reach as 
high as £p 0ak ~ 15 MeV (|Axelsson et al.ll2012D . A con- 
sideration of the usual internal shock prompt emission 
spectrum model shows that such high values are unlikely 
in this model. Furthermore, we confirm the conclusion 
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